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Abstract 

The purpose of this article is to show that the Navier-Stokes equations can be rewritten as a set of linearized 
inhomogeneous Euler equations (in convective form) with source terms that are exactly the same as those that would 
result from externally imposed shear stress and energy flux perturbations. These results are used to develop a 
mathematical basis for some existing and potential new jet noise models by appropriately choosing the ‘base flow’ 
about which the linearization is carried out. 

1. Introduction 

Lighthill (1952, 1954) provided a systematic basis for predicting jet noise when he rearranged the Navier- 
Stokes equations into the form of a linear wave equation for a medium at rest with a quadrupole-type source term, 
which includes a pressure/density contribution that Lilley (1974) showed to be more appropriately described by a 
dipole-type source. The crucial step in this so-called acoustic analogy approach amounts to assuming that the source 
term is in some sense known or that it can at least be modeled in an approximate fashion. While Lighthill’s strategy 
turned out to be remarkably successful in predicting the gross features of the sound radiation from turbulent air jets, 
the engine manufacturers soon realized that they needed a much more sensitive tool with the capability of predicting 
how even relatively small changes in the flow would affect the radiated sound. This resulted in a number of attempts 
to improve the Lighthill approach. Early efforts were focused on accounting for mean flow interaction effects. 
Phillips (1960), Lilley (1974), and many others, sought to accomplish this by rearranging the Navier-Stokes 
equations into the form of an inhomogeneous convective or moving medium wave equation rather than the 
inhomogeneous stationary medium wave equation originally proposed by Lighthill. 

But, jet noise reduction technology is still largely empirical and the most advanced industrial noise prediction 
methods, such as GE’s MGB approach (Balsa et al. 1978), which is based on Lilley’s extended acoustic analogy, 
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may do a reasonable job of predicting the noise from simple round jets, but appear to be incapable of predicting the 
changes in the sound field that occur when noise suppression devices are deployed and therefore cannot be used to 
accomplish their intended purpose — the evaluation of the acoustic performance of these devices. Tam and Auriault 
(1999) recently developed a semi-empirical model that has successfully predicted the sound at angles near 90 
degrees in a number of cases, but it does not predict the peak noise and has not as yet been used as a design tool. It 
is therefore still necessary to develop noise prediction methods that exhibit the required sensitivity to relatively 
small flow changes. 

The most obvious approach would be to use large-scale numerical simulation, i.e., DNS. However, as has been 
frequently pointed out, the number of mesh points needed to fully resolve any turbulent flow is proportional to the 
Reynolds number R based on the characteristic length scale of the turbulent eddies raised to the nine-fourths power. 
But since R is of the order of 10 5 to 10 7 for typical jet engine exhaust flows, this means that 10 12 to 1 0 15 grid points 
would be needed to resolve all of the relevant length scales. This is beyond the capabilities of currently available 
machines and the community has been forced to acknowledge that meaningful predictions will have to be based on 
some sort of combined analytical-modeling method (such as the extended acoustic analogy approaches referred to 
above) or some approximate numerical method such as LES or perhaps some combination of the two. A number of 
such approaches have recently appeared in the literature. 

The present paper is an attempt to provide a rational basis for these approaches. It begins by showing that the 
Navier-Stokes equations can always be recast into the convective form of the linearized Navier-Stokes (LNS) 
equations but with modified dependent variables, with the viscous stress perturbation replaced by a certain 
generalized Reynolds stress and, with the heat flux perturbation replaced by a generalized stagnation enthalpy flux. 
The 'base flow’ about which the equations are linearized can be any solution to a very general class of 
inhomogeneous Navier-Stokes equations with arbitrarily specified source strengths. These exact results are then 
used, along with suitable choices for the base flow source strengths, to show how many of the recent and some of 
the more frequently used classical approaches to the jet noise problem can be put on a more rational basis and, more 
importantly, extended in new directions. The LNS equations are, of course, still nonlinear but the nonlinearity is 
effectively contained in the generalized Reynolds stresses and enthalpy flux — which also contain contributions from 
the base flow sources. The focus is on obtaining useful exact results and discussion of modeling and/or numerical 
issues (and the obligatory experimental comparisons) is only included when it relates to base flow selection. 
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The acoustic analogy-type approaches and their extensions roughly correspond to treating the generalized 
stresses and enthalpy flux as known source terms that can be estimated or modeled as in the original Lighthill 
analysis. However, there are also cases where the base flow sources, which can be determined as part of the base 
flow computation, provide the dominant contribution to these quantities. This is most significant when the base flow 
is non-radiating (but, of course, unsteady) because the residual flow can then be identified with the acoustic part of 
the motion. The most obvious way of choosing the (unsteady) non-radiating base flow is to require that it be 
incompressible — but there are better choices, which are described in § 3.4 below. This leads to a different approach 
to jet nose modeling (related to the choice of base flow), which has not, as yet, been exploited. There are also cases 
where a radiating base flow is the most appropriate choice. These are described in § 3.3. 

The fundamental equations are derived in § 2 and used to justify and generalize various noise prediction methods 
in § 3. A version of the classical parallel mean flow approach (i.e. the Lilley’s equation approach) is recovered in 
section 3.1 and certain conceptual difficulties are pointed out. Section 3.2 shows how these can be overcome by 
introducing non-parallel mean flow affects. The resulting modifications can significantly change the predictions of 
the theory. 

Some of the material in this paper was presented in preliminary form in Goldstein (1999, 2000, 2002). But the 
present work is based on a more intuitively appealing form of the fundamental equations and a number of new 
results are derived. 

2. The LNS equations 

To obtain the fundamental equations, we begin with the Navier-Stokes equations in conservation law form: 
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h = h H — v 2 
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denotes the stagnation enthalpy, h denotes the enthalpy, t denotes the time, x = {.Xj,x„x 3 } are Cartesian coordinates, p 
denotes the pressure, p denotes the density, v = {v,,v„v 3 } is the fluid velocity, cr. is the viscous stress tensor, and q i is 
the heat flux vector. The dependent variables are required to satisfy an equation of state which, for definiteness, is 
taken to be the ideal gas law: 


p = pRT, h = c p T, (2.3) 

where R = c p - c v is the gas constant, c p and c v are the specific heats at constant pressure and volume and T is the 
absolute temperature. 

The dependent variables 


p = p + p', p = p + p', h = h+h', Vj — Vj + v'i , (2.4) 

as well as the viscous terms cr. and q., are divided into their ‘base flow’ components p, p, A, , and q f ,and into 

their ‘residual’ components //, //, /;', v- , <J'j , and q[ . The former are required to satisfy the inhomogeneous 
Navier-Stokes equations 
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where the nonlinear operator D is defined by 


(i6) 

for any function /and 

K=h+^v 2 (2.7) 

is the base flow stagnation enthalpy. The base flow variables, p, p, and h are also required to satisfy an ideal gas 
law equation of state, 
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and the ‘sources strengths’ Ty , H 0 , and H t , which are assumed to be localized within the jet, will be specified 

subsequently. The reason for using both over bars and tildes when defining the base flow components will become 
clear when equation (3.9) is introduced below. 

To obtain the governing equations for the residual variables, substitute (2.4) into the Navier-Stokes equations 
(2.1) and subtract out (2.5) to obtain 
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where D o (defined by (2.6)) is now a linear operator and we have put 


n = pc p T' + p'h= — p', (2.10) 

R 

and 

e ij = — (/ 7v /v) + Ty - a'lj ), (2.11) 
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Equations (2.9 b and c) can be simplified by differentiating the quadratic terms on the left hand sides by parts and 
using equations (2.5 a and b), (2.9 a), (2.10), and (2.13 b) (below) to obtain (Goldstein 1999, 2000, 2002) 
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These equations can be put into a more transparent form by introducing the new dependent variables 


and source strengths 


p'e = p' + ~~ p^’’ 2 +{r~\)H 0 , 


«i = — = P~ 
P P 


V - 1 


4 ' = e ij +L -^~ S ijP v ' 2 = -pv'p'j + L -^ L S ij pv ,A +a' ij . 


, 


,2 


and 


’I) =-pv'iK -ql + cJjjV'j, 


4 = 4 ~ S ij(r- l )Ho, 


m = H i~ T ijVj, 


to obtain 


-Dp' d - , 

P H pit i = 0, 

Dtp dxj J 


p 


f D 8vj 

U - + ll • 


Dt 1 } dx 


V 


d_ , _P^__8_~ -JLt , ~ \ 

v Pe -7. ^ v rwA e v vr 


J ) 


dxj p dx j J dx 


and 


(2.16) 

(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 

(2.22a) 

(2.22b) 


7 



y - 1 1 


D 

Dt 


Pe 


dvj 

-yp'e^+r . 

OX j OX j 


d - , 

puj 


^ T ij 8 / ~ \ / , ~ \ 5vj 


St ; St j 


(2.22c) 


where 


D _ 8 „ 8 

Dt dt ‘ dxj 

denotes the convective derivative in terms of the base flow velocity and y = c/c is the specific heat ratio. 

Equations (2.22) are the fundamental equations on which this paper is based. They are exact consequences of 
the original Navier-Stokes equations, but have been rearranged so that their left-hand sides are the same as the 
equations that would have been obtained by linearizing the convective form of the Euler equations about the base 
flow (2.5). In other words, they are the linearized inhomogeneous Euler (LIE) equations but with different 
(nonlinear) dependent variables u- and p' e , which should cause no particular difficulty since, for the present 
purpose, they can be treated as five linear equations in five unknowns — one of which reduces to the pressure 
fluctuations in the far field. The dependent variable p' e / (y - 1) is closely related to the total residual energy density 

per unit volume. The right-hand sides correspond to the sources that would be obtained by imposing the external 
stress perturbation e’j - and external energy flux perturbation - q i . In other words, the fundamental 

equations (2.22) are just the linearized Navier-Stokes (LNS) equations with modified dependent variables and with 
the viscous stress perturbation replaced by the generalized Reynolds stress ejj - e \j and the heat flux perturbation 

replaced by the generalized stagnation enthalpy flux if - fj. . A formal vector Greens function solution is given in 
the Appendix. 

Lighthill (1952, 1954) pointed out that his equation implies that there is an exact analogy between the density 
fluctuations in any real flow and those produced by a quadrupole source in an ideal (non-moving) acoustic medium. 
The LIE (or LNS) equations (2.22) show that there is an exact analogy between the fluctuations occurring in any 
real flow and the linear inviscid fluctuations about an arbitrary flow produced by appropriate external stress and 
energy flux perturbations. The externally applied stress field can be associated with a quadrupole source in the 



stationary medium case considered by Lighthill. Unfortunately, there is no completely satisfactory definition of the 



multipole source order for a general moving medium, but one possible way of generalizing the stationary medium 
definition is to associate the term quadrupole with the source produced by a fluctuating stress field. The enthalpy 
flux term would then correspond to a dipole source. This is discussed in more detail in § 3.1 below. 

The classical acoustic analogy approaches referred to in the introduction follow the Lighthill philosophy and 
rearrange the Navier-Stokes equations into the form of a single moving medium wave-like equation with an 
appropriate source term. Unfortunately, the wave operator part of the resulting equation (usually the left-hand side) 
then contains both linear and nonlinear terms. Application of these equations to specific flows requires that this 
operator be linearized by moving the nonlinear terms to the right side, which causes the source term to become very 
complicated (see Colonius, Lele, and Moin 1997 and Freund and Fleischman 2002). The present approach differs 
from these classical approaches by carrying out an exact linearization at the outset and, as the LNS (or LIE) 
equations show, the resulting source terms are then quite simple and easy to interpret — provided the right set of 
dependent variables is selected. No additional rearrangement of the LIE equations is required when specific choices 
of the base flow are introduced. 

One possible drawback of this approach is that it does not usually lead to a single wave equation (but see § 3.1 
below) for an appropriate dependent variable. This is not, however, really necessary and, as was first pointed out by 
Bechara et al. (1994), the primitive form of the equations (such as those being considered here) may even be more 
convenient for certain numerical purposes. 

3. Applications of fundamental equations 

In spite of these features, the LNS (or LIE) equations would still amount to little more than a sterile abstraction 
if specific base flow choices were not introduced. Some appropriate choices are discussed below. 

3.1 Homogeneous base flows 

We begin with the case where the base flow equations are homogeneous and inviscid so that the right sides of 
(2.5) are identically zero. A general class of solutions to these equations, which quite conveniently provides a good 
approximation to the actual mean flow field in a jet, is the unidirectional transversely sheared mean flow 

Vj = S n U(x 2 ,* 3 ), p = constant, p = p[x 2 ,x 3 ) (3.1) 
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The continuity equation (2.22 a) then decouples from the remaining four LNS (or LIE) equations, which now 


become 
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These equations are just the inhomogeneous compressible Rayleigh equations (Betchov and Criminale 1967) and, 
as is now well known (see chapter 1 of Goldstein 1976), the velocity variable u{ can be eliminated between them 
(by taking the divergence of the first equation and the convective derivative of the second, subtracting the results 
and then using the first equation to eliminate the velocity fluctuation on the left-hand side) to obtain the 
inhomogeneous Pridmore-Brown equation 
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is the variable density Pridmore-Brown (1958) operator and 


(3.6) 
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Yp! p{x 2, * 3 ) 


(3.7) 


C 2 = 


is the square of the mean flow sound speed. This is yet another form of Lilley’s equation with a quadrupole/dipole 
type source term, which is clearly more complicated than the original Rayleigh equation sources. 

This equation does not reduce to Lighthill’s result in the limit U—> 0 and P — » constant considered by Lighthill 
(1952) but it does, with some rearrangement, reduce to a form of Lighthill’ s equation obtained by Lilley (1974). 

Lighthill originally suggested that the inconvenient Sy(p — c^pj term that appears in his stress tensor should be 

associated with non-isentropic density fluctuations, which are presumably produced by viscous effects and are, 
therefore, likely to be negligible. Lilley (1974) pointed out that this term also has an isentropic component (which 
may be an important noise source when the mean density (i.e., temperature) is non-constant) and proceeded to 
derive an equation, which separated out these two components. Unfortunately, the resulting source term contains a 
dipole component (corresponding to the enthalpy flux >)'j in equation (3.3)). However, the second time derivative 
that appears in (3.6) causes it to have the same efficiency as the quadrupole source term in the low-frequency (or 
long wavelength) limit — , which is the only situation where the multipole nature of the source has any significance. 

It is worth noting that the dependent variables v/ and h' and consequently the ‘source strengths’ e-j and 77 J- 

that appear in (3.2) and (3.3) do not have zero time averages and therefore still contain an undesirable mean flow 
component. It is also worth noting that these equations possess homogeneous solutions corresponding to spatially 
growing instability waves on the base flow (3.1) (Betchov and Criminale, 1967). They also possess homogeneous 
solutions corresponding to the continuous spectrum. The complete solution to these equations consists of a particular 
solution plus these homogeneous contributions, but the result is meaningless because the instability waves become 
unbounded (infinite) far downstream in the flow. The usual resolution is to completely neglect the instability wave 
contribution. Dowling, Ffowcs Williams, and Goldstein (1978) attempted to justify this by invoking causality 
arguments. 

3.2 Steady non-parallel base flow 

A better resolution is obtained by choosing the base flow to be the actual mean flow of the jet. The over bar on 
the dependent base flow variables would then denote the time average 
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denotes a Favre averaged quantity (Lele 1994) for all variables except h Q , which is defined by (2.7). Notice that 
(2.8) is completely consistent with the overall ideal gas law (2.3) when the tilde is defined in this fashion. 

The time derivatives drop out of the base flow equations (2.5), which do not, of course, form a closed system. 
The source strengths are given by 
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The base flow equations are now the ordinary RANS equations, which are usually closed by assuming some 


sort of model relating the source terms to the mean flow variables vj- , P , P , and h and their derivatives, such as 


the Boussinesq model (Speziale 1991; Speziale and So 1998) 
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for the Reynolds stresses, with additional equations to determine the turbulent viscosity p T and a similar model for 


Hj . The fundamental LNS equations (2.22) remain exact even when these approximations are introduced since, as 
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already noted, the base flow source strengths can be arbitrarily specified. But, (3.8) would only be an approximation 
and the base flow would be somewhat different from the actual mean flow in the jet. 

The time average flow need not, of course, be determined from a RANS computation. It can, for example, be 
determined by applying (3.8) to a LES or a DNS computation. The base flow sources (3.10) and (3.12) can then be 
written more compactly as 


fj=-pv\v , j , (3.14) 

V i =H i -T ij V j =-p (3-15) 

with H 0 still be given by (3.1 1) and ^ h 0 -h$ is given by (2.15). The inviscid LIE (or LNS) equation source 
strengths e'y - % and rj ■ - //, would be proportional to the differences between the fluctuating and Favre averaged 

Reynolds stresses and enthalpy fluxes, which means that these sources would now be true fluctuating quantities 
with zero time averages, even in the viscous case .This insures that the residual variables will have that property as 
well. (The fact that p' e =u' = 0 can also be verified directly from (2.16), (2.17), (3.11), and (3.14).) 

The corresponding LIE equation sources, which are in principle exact and (as noted in section 2) simpler than 
those that appear in the literature, must still be modeled. This is most easily accomplished by using the formal 
Greens function solution (A-9) to express the far field pressure in terms of the near field source distribution. The 
result can then be time averaged to relate the mean square pressure in the far field to the source correlation function, 
which means that it is only necessary to model this latter quantity and not the instantaneous sources themselves. The 
resulting fourth order two point time delayed correlations involve pure fluctuating quantities and do not contain the 
awkward mean flow component that appears in the parallel base flow model discussed in § 3.1. (The importance of 
this property was demonstrated by Bogey, Bailly, and Juve (2002)). These correlations can be expressed in terms of 
second order velocity and enthalpy correlations (which are very close to what is experimentally measured) by 
invoking the quasi-normal approximation (Millionschikov, 1941) — which is thought to be fairly accurate in this 
context. 

The most important difference between the present results and the parallel flow results of § 3.1 is that the 
homogeneous solutions to the LIE (or LNS) equations, which presumably correspond to instability waves growing 
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and then decaying on the diverging non-parallel base flow, will always remain bounded. This eliminates the paradox 
that arose in § 3.1. The complete LIE equation solution therefore remains bounded and can, as a consequence, be 
used to calculate the radiated sound. It consists of a particular solution that is driven by the sources (i.e. it satisfies 
causality and therefore the appropriate upstream boundary conditions) and a homogeneous solution that is driven by 
the upstream boundary conditions, say at the nozzle exit. 

The significance of this result can be understood by considering an idealized subsonic or a nearly perfectly 
expanded supersonic jet. The spread rate, say e, of the base flow will then be uniformly small, which implies that the 
mean pressure gradient will be 0(e), the mean cross-flow will be 0(f), and all streamwise changes in the base flow 
will occur on the slow streamwise length scale sx v The source driven particular solution can, as noted above, be 
constructed from the casual vector Green function, which satisfies (A-7) and consists of a component that decays to 
zero when the unsealed streamwise coordinate x x - .%• becomes large plus a component that becomes unbounded 
when the unsealed streamwise coordinate becomes large but decays to zero on the long (slow) stream-wise length 
scale £ ( x, - Xj*). The former component can be calculated by treating the slow variable 8 x , as a parameter and 
using the parallel flow solution (3.1) to simplify the linearized Euler operator (A— 1). 

This approximation cannot, however, be used to determine the latter component (which corresponds to a 
superposition of linear instability waves, Briggs, 1964; Bers, 1975), since the resulting local solution would become 
unbounded and would not remain valid over the long streamwise length scale £ ( jc, - x*) on which this component 
evolves. The appropriate result is obtained by retaining 0( £ ) terms in the linearized Euler operator (A-l) and using 
the multiple scales method together with the WKBJ approximation to render the solution uniformly valid. 

Equation (A-8) (or (A-9) and (A- 10)) shows that both components of the Greens function act on the same 
source term. Each component is capable of producing acoustic radiation - even at subsonic Mach numbers (Tam 
and Morris 1980). The first corresponds to the usual Lilley equation solution but with slowly varying coefficients 
and with e'y and >l'i replaced by e\j - % and ifi ~ Hi , respectively. The second is associated with linear instability 

waves but is very different from conventional instability models, since these waves are now continuously generated 
along the length of the jet and do not constitute separate sound sources. Their only role is to provide another path 
for propagating the source generated pressure fluctuations to the far field. 

Each Greens function component can be thought of as a filter that only allows certain elements of the source 
spectrum to reach the far field. Since these components are expected to respond to different portions of this 
spectrum, the resulting radiation field should exhibit a bi-model structure. Tam, Golebiowski, and Seiner (1996) 
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recently presented some experimental evidence which strongly suggests that the jet sound field actually has such a 
structure — especially at supersonic Mach numbers. The contribution of the second Greens function component will 
probably be fairly small at subsonic Mach numbers but is likely to be dominant at supersonic speeds. The present 
result is therefore very different than the parallel flow model of the previous subsection. 

The homogeneous solution, which is usually associated with the large-scale coherent structures (Crow and 
Champagne 1971; Mankbadi and Liu 1984; Gaster, Kit, and Wygananski 1985), is similar to the second component 
of the Greens function solution in that it involves instability waves and must also be determined from a multiple 
scales analysis. It is, however, expected to be small compared to the contribution of that component, since it is only 
generated at the upstream boundary while the latter is generated over an extended streamwise region. However, the 
only unequivocal distinction between the two solutions is that the homogeneous solution generates sound through a 
linear process while the particular solution generates sound through non-linear mechanisms. They can, of course, be 
coupled since the acoustic radiation from the particular solution can trigger instability waves in the upstream region, 
say at the nozzle lip. But, this effect is probably insignificant (relative to the effect of the pre-existing upstream 
disturbances) for the noisy flows that are of technological interest. 

Another difference between the present results and those of the previous subsection is the r ( y term on the left- 

hand side of (2.22 b), which drops out of the parallel base flow model of § 3.1 but can now produce a potentially 
strong coupling between (2.22 b) and (2.22c) and the continuity equation (2.22a) . This term is expected to be small 
when the mean flow is nearly parallel, as it is in idealized subsonic jets. However, real turbojet flows are not all that 
parallel in the vicinity of the nozzle — especially when the latter is of the mixer-ejector type. The term could also be 
important for improperly expanded supersonic jets and in flows that arise in non-aeronautical applications of the 
LNS equations (2.22). 

3.3 Unsteady compressible base flows 

While the coherent structures are at best poorly defined, there is some evidence to suggest that they may not be 
adequately represented by the linear instability waves. A better representation can be obtained by making them part 
of the base flow, which would then, of course, have to be unsteady. Calculations based on linear stability theory 
(Reau and Tumin 2000a and b; Liu 1998) suggest that good results can be obtained by allowing the base flow 

momentum source Ty to again be of the Boussinesq form (3.13) with a similar representation for the enthalpy 
source Hj . H a would still be given by (3.18). The RANS type eddy viscosity model also seems to work fairly 
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well. The resulting base flow equations would then be the unsteady RANS or (URANS) equations with the 
unsteadiness being driven by the random upstream boundary conditions. The LNS (or LES) equations (2.22) are 
again exact. 

Bastin et al. (1997) carried out the relevant computation for a two dimensional planar jet. They used a standard 
k-s model for the eddy viscosity with an appropriately reduced numerical coefficient. Their results seem to over 
predict the contributions of the large scale structures, which suggests that a more refined turbulence model might be 
in order. 

The relation between the URANS solutions and the coherent structures can be made more plausible by letting 
the over bars on the base flow variables denote the ensemble/phase average 

M N 

f = lim V V fj (x,t +/,-.) (3.16) 

oo NM ^ ^ 1 v tJ ’ v ’ 

7=1 t=l 

where f can be any of the dependent variables, f denotes the value of f in the i ,h realization of the flow and the t.. 
denote appropriate time delays, which are introduced to account for the fact that each of the M coherent structures 
will have a different phase in each realization of the flow. Bastin et al. (1997) proposed using a conventional 
ensemble average to represent the coherent structures, which would correspond to setting r = 0 in. (3.16). Equation 
(3.16) will represent a conventional phase average when M=1 and all of the members of the ensemble correspond to 
the same realization of the flow. This would probably provide a good representation of the coherent structures in a 
flow that is artificially excited by a periodic source (Liu, 1998). The appropriate time intervals t u would then be 
integral multiples of the period T of the excitation. The precise values of the r would, of course, be unknown in 

naturally occurring flows, but that is largely irrelevant for the present purposes. / consists of the time average of f 
(i.e., the mean flow) plus the ensemble/phase average of the fluctuating component of / (Reynolds and Hussain 
1970; Hussain 1983). 

The tilde is still defined by equations (2.8) and (3.9), the base flow source terms Tj , H 0 , and If are still given 

by (3.10) to (3.12) (but not by (3.14) and (3.15)) and the base flow is no longer non-radiating. It is likely to be 
composed of the larger scale components of the motion and may therefore bear some resemblance to the resolved 
component of a corresponding Large Eddy Simulation (LES) (Bastin et al. 1997). As in LES, the sources, which 
now account for the effect of the random turbulence on the coherent motion, must be modeled. Using the 
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Boussinesq type model (3.13) for the momentum source (3.10) and a similar model for the enthalpy source (3.12) 
leads to the URANS equations alluded to above. 

While still expensive, URANS computations are much less demanding of computer resources than the DNS 
computations referred to in the INTRODUCTION (Bastin et al. 1997). There is, therefore, some hope that they can 
be used to predict the sound emission from realistic flows. In any event, the residual variables (which now 
correspond to the incoherent part of the motion) can, and probably do, contain significant large scale components, 
which would not be the case in LES where the pure spatial filtering that is usually used would cut off these 
components. These variables are again governed by the LIE (or LNS) equations (2.22) but the homogeneous 
solution (if it exists) must now be set to zero, since the upstream boundary conditions that would drive this solution 
are already satisfied by the base flow. This approach formally resembles the two stage LES approach proposed by 
Bodony and Lele (2002), but the philosophical basis and modeling assumptions are quite different. 

The complete sound field can, at least in principle, be calculated by adding the base flow sound field to the 
far-field LIE (or LNS) equation solution. The latter equations are, in principle, exact, but their source terms still 
have to be modeled .It is, however, no longer possible to use the procedure described in the previous subsection to 
relate the far-field pressure correlation to the correlation of the source function since the base flow is now unsteady. 
It is therefore necessary to model the instantaneous source strengths rather than their correlations — which is a much 
more difficult task. It may be possible to use some sort of stochastic source model such as the ones used by Bailly et 
al. (1995) and Bodony and Lele (2002). A more expedient approach might be to neglect the contribution of the 
unsteady coherent component of the base flow to the coefficients of the LIE equation operator ( A— 1), since it is 
likely to be fairly small compared to the mean flow contribution. The incoherent far-field pressure auto-correlation 
could then be related to the source correlations, which would now, however, involve only the incoherent part of the 
motion. But this creates somewhat of a dilemma because only the correlation of the total unsteady motion is actually 
measured in the experiments. And even if this difficultly could be overcome, the result would only be meaningful if 
the coherent and incoherent components of the motion were relatively uncorrelated (Bastin et al 1997). It would then 
be a generalization of the RANS equation approach described in the previous subsection and would reduce to that 
approach if the base flow computation were linearized about the mean flow and source term perturbations were 
neglected. 

However, we have already indicated that this linearization may be unrealistic and Reau and Tumin (2002a and 
b) show that the perturbations in the source terms cannot be neglected. Liu (1988) and Mankbadi and Liu (1984) 
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adopted an intermediate approach in which the coherent motion was determined from integral conservation laws by 


using an approximate shape functions taken from linear instability theory. 

The base flow equations with source strengths (3.10) to (3.12) formally resemble the filtered Navier-Stokes 
equations (Rogallo and Moin, 1984) and there are a number of interesting ways of choosing the filter in these 
equations. Goldstein (2002) showed that it can even be chosen to only pass wave number components with subsonic 
phase velocity in the stream-wise direction* which, according to conventional wisdom (see, for example, Freund 
2001), should cause the unsteady base flow to be non-radiating .But this probably needs further study. The relevant 
flow can be computed by using a Fourier-Spectral method with only subsonically converting spectral base 
functions retained in the computation. The advantage of using an unsteady but non-radiating base flow is discussed 
in the next subsection. 

3.4 Non-radiating unsteady base flows 

The radiated sound is typically many orders of magnitude smaller than the non-radiating component of the 
motion in virtually all high speed air jets, as well as in most other high speed flows. In fact, it is usually many orders 
of magnitude smaller than the errors incurred in computing the non-radiating part of the flow. A number of 
investigators have even suggested that the ‘acoustic component’ of the numerical solution to the full Navier-Stokes 
equations would be hopelessly corrupted by these errors or even by the computational noise itself (e.g. Roe 1992; 
Crighton 1993). This difficulty would certainly be overcome if the radiating component of the motion could be 
separated out and calculated as if it were a small perturbation about the non-radiating component, which amounts to 
requiring that the base flow in the decomposition (2.4) be the entire (or almost the entire) non-radiating component 


*There is a subscript 1 missing from the k on the right side of equation [49] of that reference, 
of the motion. 
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A nice feature of the RANS base flow of § 3.2, or for that matter, any steady base flow, is that it is non- 
radiating. Another possible non-radiating base flow is an incompressible flow. More generally, the base flow can 
always be made non-radiating by requiring that 

-^ = o(|x| 4 j as |x_l | — » oo, (3-17) 

where jc ± is the transverse coordinate vector. This will force the flow to behave incompressibly at large transverse 
distances and should therefore prevent any acoustic radiation from reaching the far-field. 

It is, in general, impossible to satisfy the homogenous energy equation (2.5c) and still impose arbitrary 

requirements on op / of ■ We can, however, choose the energy equation sources H 0 and 1I 1 so that this equation is 
identically satisfied for any p by setting 

(y-l)^o =Pe =P + -^~Pv 2 , =pv i +q i -o j v j (3.18) 

The momentum equation source strength 7 )j can again be set to zero. It is then possible to preselect p and require 

that is satisfy (3.17). This means that (2.5a) and (2.5 b) are similar to the incompressible Navier-Stokes equations in 
that they can be solved for vj- and p independently of the energy equation. 

The source terms in the LNS equations (2.22 b) and (2.22c) now become 


e ij e ij “ ' 


V — 1 9 | f _ y — l 9 

-^ s ii v ^Y a [] +s ii\v+—^py 


vqt'j 


(3.19b) 


Vi ~ Vi = -pv'iK - pvf - q t + GjVj - cr'jjVj 


(3.19c) 


A cleaner result is obtained when H 0 is eliminated from the residual equation dependent variable p' e by letting 


Pe ^Pe -(r-l)H 0 =P' +^—^-pv' 2 - 


(3.20) 
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The left sides of the LNS (or LIE) equations (2.22) remain unchanged and their right-hand sides become 


RHS (2.22a) = 0. 

(3.21a) 

dc'- 

RHS (2.22b) = — , 

V ’ dx : 

(3.21b) 


and 


RHS (2.22c) -> - 


f D 




y - 1 


— V + YV 

Dt 


5x 


J J 




Ox ; 


(3.21c) 


The strictly incompressible base flow is recovered by setting p =constant. Hardin and Pope (1994) developed a 
two-stage theory in which the base flow was incompressible and viscous but the residual flow was homentropic. A 
more consistent approximation would be to set the residual viscous and heat conduction terms to zero in (3.21 b and 
c), which roughly corresponds to the somewhat weaker requirement that the residual flow be isentropic. 

When the flow is inviscid and the base flow density p is constant, the deviation v\ of the velocity form and the 
incompressible approximation vj- should be 0 (M ) when the Mach number M is small and the quadratic source terms 

d{e\j / ox j j and [drj'j / dxj j + e’y (dv i / dx i ) should then be negligible compared to the base flow source term in 

(3.21 c). This implies that the sound field is primarily generated by the convective derivative of the incompressible 
pressure fluctuations, since dv i /dx i = 0 in this case. The result is, therefore, related to Ribner’s (1962) ‘dilitation’ 

theory of jet noise, but unlike Ribner’s original proposal, the interaction of the sound with the incompressible flow is 
now explicitly accounted for. Also, the purely convected pressure disturbances (i.e., the pseudo sound) are now 
(explicitly) non-radiating. The important thing is that the sound sources are nearly independent of the residual 
variables, which should certainly be required if the latter are to be associated with the acoustic part of the motion. 
Tam and Auriault (1999) have also proposed using the pressure fluctuation (though not the incompressible pressure 
fluctuations) as a sound source. 

The strictly incompressible approximation cannot remain valid at supersonic Mach numbers (even for ideally 
expanded jets), since shear layer spread rates and turbulence structures are known to be Mach number dependent in 
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this case (Smits and Dussauge 1996, ch. 6), and it is even questionable at moderately high, say above 0.3, subsonic 


Mach numbers. It may, however, be possible to choose the density p so that the base flow (2.5a) and (2.5 b) with 
Tjj = 0 remains close to the real flow and it still non-radiating (i.e., it satisfies (3.17)). An intuitively appealing way 

of accomplishing this is to require that p be convected by a localized flow field, say V t ( x , / ) , whose streamlines are 
constrained to lie within the interior of the jet, which means that 


d_ 

dt 



+ V; 


3x : 



= 0 , 


(3.22) 


where R becomes constant at large distances from the flow, but is otherwise arbitrary. The upstream boundary 
conditions can then be chosen to insure that P decays at large transverse distances from the flow, which should then 
ensure that (3.17) is satisfied. This can be thought of as a generalization of Morkovin’s hypothesis (Morkovin 1962; 
Bradshaw 1978) which, unlike the original proposal allows for fluctuations in the density. There is therefore some 
hope that this model (unlike the original Morkovin proposal) will provide a reasonable approximation for jet-like 
flows, even at supersonic Mach numbers (Smits and Dussauge 1996, p. 138). 
p and Vj will also satisfy an equation of the form 




(3.23) 


if Vj and R are required to satisfy the steady flow type continuity equation 


— RV- =0 
5xj 


(3.24) 


In which case, the continuity equation (2.5a) becomes 
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(3.25) 


a 

SX; 



Vi) = 0, 


which means that the base flow satisfies the steady continuity equation in a reference frame moving with the 
velocity Vj , i.e., the relative momentum flux p ( '■’/ - V i j now has zero divergence. 

For numerical purposes it may be better to determine p from an equation of the form 


dp Sjj 

at a Xi a X j 


(3.26) 


where Sy (x, t) is required to vanish at an appropriate rate as .r— »oo, but is otherwise arbitrary. This can be seen by 

eliminating the momentum flux ~pv t between the base flow continuity and momentum equations (2.5a) and (2.5 b) 
to obtain the Poisson’s equation 

7 _ a 2 s’;* a^ViVj-^+s.j) 

p = = 

dx^dxj dx^dxj 



where the dot denotes differentiation with respect to time. The free-space Greens function solution to this equation is 


p(x.t) 


4 K 


a 2 P^'') 

dxjdx j J \x' - x\ 


(3.28) 


which shows that p = o(|x| 3 j as |x| — > oo if Sy is chosen so that Sy = o||x| 3 j as |x| — > oo . This is the same as the 

pressure decay rate experienced by a localized incompressible flow (Crow 1970) and is therefore likely to be fast 
enough to greatly simplify the implementation of the far-field boundary conditions. The density p will still satisfy 
(3.23) (and therefore (3.25)) if is defined by 
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(3.29) 


V ; 


1 5§ ij 

P 5x J ’ 


The density can easily be localized if V t is required to be steady. It could even be taken to be the actual mean 
velocity of the jet, as determined from a preliminary RANS calculation (see § 3.2). In which case, (3.24) would 
require that R be the mean density of the jet. 

There are, of course, many other possible choices for R and V: ■ They could, for example, be chosen so that the 
base flow source term 


f TT 


V - 1 


D _ -dvj 

a p * rp & T 


J ) 


in (3.21c) has zero mean, which would then ensure that the base flow satisfies the correct homogeneous energy 
equation in the mean and that p' e has zero mean. This would, however, require that an additional integral- 
differential equation be solved along with (2.5a and b), which would make the base flow computation more 
difficult — possibly more difficult than an actual compressible DNS computation, which would, of course, defeat the 
whole purpose of the two stage approach. 

Lighthill (1952) argued that the strength of his quadrupole source could be obtained to a good approximation by 
calculating its value for an equivalent flow devoid of sound. The present result provides an analytical basis for that 
idea. It is, of course, possible to move the residual stresses to the left side of the LNS equations and calculate the 
sound from the full nonlinear equations, but that would make the present approach more complicated and 
computationally more expensive than solving the original Navier-Stokes equations. In fact, the main justification for 
using this approach is to ensure that the errors incurred in computing the non-radiating component of the motion do 
not unduly contaminate the radiating component — which is best accomplished can only be accomplished if the 
residual variables are required to be much smaller than the base flow variables, which should ensure that the 
linearization is justified. 


23 



4. Conclusions 


A general set of linearized inhomogeneous Euler (LIE) equations has been derived and used to put a number of 
classical and recently proposed noise prediction methods on a more rational basis and in most cases extend them in 
new directions. The source terms are exactly the same as those that would be obtained by imposing externally 
applied stress and energy flux perturbations on the flow and the equations are therefore exactly the same as the 
linearized Navier-Stokes (LNS) equations, but with the viscous stress perturbation replaced by an appropriate 
Reynolds stress and the heat flux perturbations replaced by an appropriate enthalpy flux. The dependent variables 
are also somewhat different. 

Lighthill pointed out that his equation implied that there was an exact analogy between the density fluctuations 
in any real flow and those produced by a quadmpole source in an ideal (non-moving) acoustic medium. The LNS (or 
LIE) equations show that there is an exact analogy between the fluctuations in any real flow and the linear inviscid 
fluctuations about an arbitrary base flow produced by externally imposed stress and energy flux perturbations. This 
does not, of course, imply that the LNS (or LIE) equations, or for that matter, any of the previously proposed 
‘acoustic analogy’ equations can provide an unambiguous identification of sources, except perhaps for the 
generalized incompressible flows considered in § 3.4. These equations are only useful when the ‘base flow’ is 
reasonably close to the actual fluid motion and, in most cases, can only serve as a guide for identifying and 
ultimately modeling the apparent sources of sound. The sound is just a by-product of all the processes occurring in 
the flow and it is highly unlikely that “true” sound sources can be identified in any realistic turbulent flow. 

The LIE equations are made explicit by considering various possible base flows. A number of such flows, 
corresponding to progressively more complex choices of the base flow source terms, were described in § 3. 

The first three choices require source term modeling as in the classical Lighthill approach. But these terms now 
generate only part of the sound field in at least two of these. The remainder must be determined by the direct 
solution of the appropriate linear or non-linear equations, which is another indication that the mathematical sources 
need not be the true sources of sound. The last subsection considers formulations in which the dominant LNS 
equation source term can be computed but require that the base flow be chosen to minimize the neglected source 
terms. — which is a different approach to jet noise modeling that has yet to be exploited. 
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Appendix 


The five LNS equations (2.22) can be written more compactly by introducing the five dimensional operator 


L 


/uv 


8 /xv D o +(r- !)^ v 4 + d 1 


r - 1 


+^t/5 


+ K 


MV 


(a - 1) 


where //, r= 1,2,. . . ,5 while the Latin indices i,j are restricted to the range 1,2,3. 
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1 dr 


v4 ~- 


V] 


P dx 


°M 4 


J ) 


(A -2) 


8ju -= 


— , i = fi= 1,2,3 

OX; 


0, otherwise 


(A -3) 


v m and Tj are both zero when /j> 3, c 2 is given by (3.7), the Einstein summation convention is still being used 
and S Llv is the five dimensional Kronecker delta. 

Then Eq. (2.22) can be written as 


L mv u v 


(A -4) 


where 


«v =i m i’ 


Pe 

y-l : 


(A-5) 


is the five dimensional dependent variable vector, the five dimensional source vector s is given by 
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(A -6) 


dx j 


(' e Mj e M j ) + 3 /j4 ( e 'tj 



and s\j = t]'j and £4 ,■ = // ; . This system can be formally solved by using the vector Greens function (Morse and 
Feshbach 1953, pp. 878-886) g va ( x d , which satisfies 

L juv gva = - x Vi 1 - 1') ’ (A - 7) 


to obtain 





—00 


J g (*, *|x V 'df' 
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(A -8) 


and the derivatives can be transferred to the Greens function to obtain 


/ (*> 0 = - J J r vj,n (e’tf - )dx 'dt ' 


(A -9) 


-00 V 
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